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Abstract 

In his deep and prolific investigations of heat diffusion, Lame was led to the inves- 
tigation of the eigenvalues and eigenfunctions of the Laplace operator in an equilateral 
triangle. In particular he derived explicit results for the Dirichlet and Neumann cases 
using an ingenious change of variables. The relevant eigenfunctions are complicated 
infinite series in terms of his variables. 

Here we first show that boundary value problems with simple boundary conditions, 
such as the Dirichlet and the Neumann problems, can be solved in an elementary man- 
ner. In particular for these problems, the unknown Neumann and Dirichlet boundary 
values respectively, can be expressed in terms of a Fourier series. Our analysis is based 
on the so called global relation, which is an algebraic equation coupling the Dirichlet 
and the Neumann spectral values on the perimeter of the triangle. 

As Lame correctly pointed out, infinite series are inadequate for expressing the 
solution of more complicated problems, such as mixed boundary value problems. Here 
we show, utilizing further the global relation, that such problems can be solved in terms 
of generalized Fourier integrals. 



1 Introduction 



Solutions of certain linear elliptic boundary value problems can be expanded in complete 
sets of eigenf unctions. Unfortunately, the actual form of these eigenfunctions is known for 
only simple geometries. In fact, only geometries that allow separation of variables yield well 
known expressions for the associated eigenfunctions. But what happens when separation 
of variables does not apply? Is it possible to construct the spectral characteristics of a 
fundamental domain that does not fit any separable coordinate system? Some examples 
where this construction is possible are presented in the present work. The approach used 
here has its roots in the unified transform method for analysing both lineal and integrable 
nonlinear PDEs introduced in [4]. 

A crucial role in this analysis is played by a certain equation coupling all boundary 
values, which was called the global relation in [4]. The concrete form of this equation for the 
equilateral triangle was given in the important work of [14], where it was called a functional 
equation. 

A general overview of the problems solved in this paper is presented in the sequel where 
notations and some elementary formulae are included in order to facilitate the understanding 
of the new results. We study boundary value problems for the Laplace, the Helmholtz and the 
modified Helmholtz equations in the interior of an equilateral triangle. These equations are 
three of the basic equations of classical mathematical physics. In particular, they arise as the 
reduction of several fundamental parabolic and hyperbolic linear equations. Furthermore, 
the specific boundary conditions discussed here cover most cases of physical significance. 

We first introduce some notations. 

1.1 Notations and Useful Identities 

(i) z will denote the usual complex variable and a will denote one 
unity, 

2iTT 1 iy/S 

Bar will denote complex conjugation, in particular 

z = X — iy, a = e 3. 

F(k) will denote the Schwarz conjugate of the function F{k). 
(ii) The complex numbers 

I ijl _ I , . 

Zl = -j=e 3 , Z2^ Zl, zs = --j=, (1.2) 

will denote the vertices of the equilateral triangle, and D G C will denote the interior of the 
triangle. The length of each side is /. 

The sides {z2,zi), {zs,Z2), (zi^zs) will be referred to as sides (1), (2), (3) respectively. 



of the complex roots of 
(1.1) 
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Figure 1.1: The fundamental domain D 

(iii) On each side we identify the positive direction T and the outward normal N as in 
Figure 1.1. The functions 



I I 
2' 2 



J = 1,2,3, 



:i.3) 



will denote the function q{x, y), as well as the derivative of q{x, y) along the outward normal 
N respectively, for the side (j). 

(iv) E{k) and e{k) will denote the following exponential functions 



E{k) = exp <^ {k + T 



A 



(v) Using the fact that the numbers a and a satisfy the obvious relations 

= a = a~^, 1 + a + a = 0, ia — ia = \/3, ia — i = \/3a, 



[lA) 



[1.5) 



it is straightforward to obtain analogous relations for E{k) and e{k). For example, the last 
three equations in (1.5) imply 

E{k)E{ak)E{ak) = 1, E{iak)E{-iak) = e{k), E{iak)E{-ik) = e{ak). (1.6) 



1.2 Formulation of the Problem 

We will investigate the basic elliptic equations in the interior of the equilateral triangle D, 
namely we will study the equation 



qxx + qyy - 4Ag = 0, (x, y) G -D, 



:i.7) 



where q{x,y) is a real valued function and A is a real constant. The case of A = 0, of A 
negative, and of A positive, correspond to the Laplace, the Hclmholtz, and the modified 
Helmholtz equations, respectively. We will analyze the following problems: 
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(i) The Dirichlet problem 



l_ l_ 
2' 2 



J = 1,2, 3. 



(1.8) 



(ii) The obhque Robin problem, 



s\n(5q%\s) + cos (3-^q^^\s) + 'yq^^\s) = fj{s), s G 



2' 2 



J = 1,2,3, 



:i.9) 



where /3 and 7 are real constants and sin /? 7^ 0. The sum of the first two terms of the Ihs 
of this equation equals the derivative of q^\s) in the direction making an angle /3 with the 
positive direction of the side (j). The Neumann and the Robin problems correspond to the 
following particular choices of /? and 7, 



Neumann : /3 = — , 7 = 0; Robin : = — , 7 7^ 0. 



(iii) The Poincare type problem 



sin (s) + cos (3j ^q^'^ (s) + 'yjq^'\s) = (s), se 



I I 
'2' 2 



J = 1,2,3, 



:i.io) 



where Pi is a real constant such that sin /3i 7^ 0, P2 and /?3 satisfy sin /32 7^ 0, sin 7^ and 
are given in terms of /3i by the expressions 



TITT TTITT 

/92 = A + -3-, '^3 = /3i + -3-, m,neZ, 
and the real constants {7^}? satisfy the relations 

sin3/?i [72(3A - jI) - e^"^7i(3A - 7^)] = 0, 
sin3A [73(3A - 7I) - e^"^'^7i(3A - 7?)] = 0. 



;i.l2a) 

(1.126) 
(1.12c) 



A particular case of such a Poincare type problem, which is solved in detail, is the modified 
Helmholtz equation with Neumann values on sides (2) and (3) and with Robin values on 
side (1), where the constant 7 is given by VSX. 

We assume that the functions fj have sufficient smoothness and that they are compatible 
at the corners of the triangle. The case of boundary conditions which are discontinuous at 
the corners will be considered elsewhere. 



1.3 The Global Relation 

As it was mentioned earlier the approach used here is based on the analysis of the global 
relation, which is the fundamental algebraic relation that couples the Dirichlet and the 
Neumann values around the perimeter of the triangle. This equation, first derived for the 
case of equilateral triangle in [14] (see also [5]) is 

E{-ik)^i{k) + E{-iak)^2{ak) + E{-iak)^3{ak) 
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^2i{E{-ik)^i{k) + E{-iak)^2{o(k)+E{-iak)^3{ak)}, A; e C - {0}, (1.13) 

where the exponential function E{k) is defined in equation (1.4a), and and $j are the 
following transforms of the Neumann and Dirichlet boundary values: 



^j{k) 



exp 



exp 



(1.14) 

for each j — 1,2,3, and every complex k ^ 0. 

The general methodology introduced in [4], [5] implies that the global relation must be 
supplemented by its Schwarz conjugate, as well as by the four equations obtained from these 
two equations by replacing k with ak and with ak. We will refer to these six equations as 
the basic algebraic relations. In this paper we present two different techniques for solving 
these equations. 



ds, 



1.3.1 Solutions via Infinite Series 

For simple problems it is possible to compute the unknown boundary values by evaluating 
the basic algebraic relations at particular discrete values of k. This yields the unknown 
boundary values in terms of infinite series. The Dirichlet and the Neumann problems are 
examples of problems which can be solved using this technique. 

We use the Dirichlet problem to illustrate this approach. In this case the functions 
appearing in the rhs of the global relation (1.13) can be immediately computed in terms of 
the given boundary conditions fj, thus the global relation becomes a single equation for the 
three unknown functions {^jji- Multiplying this equation by E{iak), and multiplying the 
Schwarz conjugate of the global relation by E{—iak), we find the following two equations 
(where we have used the last two of the identities (1.6)) 

e{ak)<i/i{k) + e{-k)^!2{ak) + ^3(0;^;) = 2iA{k), (1.15) 

e{-ak)^!i{k) + ^2(0^) + e{k)^!z{ak) = -2iB{k). (1.16) 

In these equations A{k) and B{k) are known functions and /c e C — {0}. 

For the general Dirichlet problem, we will supplement these two equations with the four 
equations obtained from these equations by replacing k with ak and with dk. However, there 
exists a particular case for which it is sufficient to analyze only the above two equations. This 
is the symmetric Dirichlet problem, namely the problem where the functions fj are all the 

same, /-,• = /, j = 1, 2, 3. Then the Neumann values q^^\s) are also the same, q^^^ = qN and 
hence ^j{k) = "^{k), j = 1,2,3. Thus equations (1.1 5) and (1.16) become two equations for 
the three unknown functions ^{k), ^{ak), ^{ak). Hence, any two of them can be expressed 
in terms of the remaining one, for example ^(aA;) and ^(afc) can be expressed in terms of 
^(/c). In particular, subtracting equations (1.15), (1.16), we find 

{e{k) - e{-k)) ^{ak) = {e{ak) - e{-ak)) ^{k) - 2iG{k), (1.17) 

where G{k) = A{k) + B{k) is a known function. Equation (1.17) is a single equation for 
the two unknown functions 'if{ak) and ^'(A;). However, by evaluating this equation at those 
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values of k for which the coefficient of ^{ak) vanishes, i.e. at e^(A;) = 1, or A; = s„, 



A 2imi 
Sn-\ = — n e Z, 

Sn. I 



(1.18) 



it follows that ^(s„) can be determined. Recalling the definition of ^(/c) and evaluating 
equation (1.17) at k = Sn, we find 



sinh 



aSr, 



X 



aSr 



qN{s)ds = iG{Sn), n G Z. 



:i.i9) 



Thus q'Ar(s) can be expressed as a Fourier series. 

For the general Dirichlet problem (1.8), the six basic algebraic relations couple the nine 
unknown functions {^^(A;), 'ifj{ak), ^^(aA;)}^. Thus any six of them can be expressed in 
terms of the remaining three. In particular, it is shown in section 3 that ^2{oik) can be 
expressed in terms of {'^j{k)}l by the equation 

{e^{-k) - e^{k)) *2(aA;) = [e{-ak) - e'^{-k)e{ak)] {^i{k) + e^{k)^3{k)) 

+e\-k) [e{-ak) - e\k)e{ak)] ^2(A:) + 2iX{k), (1.20) 

where X{k) is known. In spite of the fact that this equation is a single equation for four 
unknown functions, it yields all the three Neumann values q^j^\ j = 1,2,3. Indeed, by 
evaluating equation (1.20) at those values of k for which the coefficient of \E'2(q;A;) vanishes, 
i.e. at e^{k) = 1, or k — k^, where 



A 2irmT 
km + — ^ ^, , m e Z, 



kn 



31 ' 



equation (1.20) yields 



e 3i - 

I 
1 



g^^(s) + e-"^g)v'(5) + e"^g]v'(s) ds = M{k^), me'L, 



,(2) 



(3), 



(1.21) 



;i.22) 



where M{km) is known. This equation in contrast to equation (1.19) involves three unknown 
functions. However, equation (1.21) gives three times as many values for m as equation 
(1.18). Replacing in equation (1.22) m by 3n, 3n — 1, 3n — 2, and inverting the left hand 
sides of the resulting equations, we find 



e 3' 



e 3! 



— oo 

oo 



q'i\s) + aqPis) + aq^^\s) = y ^ M(3n - 2)e 



,(3) 



(1.23) 



Thus by solving this system of three algebraic equations it follows that each one of the 
Neumann boundary values can be represented in terms of a Fourier series (see Proposition 
3.2). 

The analysis of the oblique Robin problem (equation (1.9)) is similar. However, the 
values of s„ and of km in general cannot be found explicitly. The values k^ satisfy the 
transcendental equation 

ei'-^ity^ ^ ^ \ , ( =e^, meZ. (1.24) 

[^kme-'^ + - 7) (akme^^ + ^ - 7) 

Thus in equation (1.22) instead of exp{^} we now have exp | km + ^ | > where km. satisfies 
(1.24). In the particular case of the Neumann problem, km satisfies equation (1.21). 



1.3.2 Solutions Via Generalized Fourier Integrals 

For more complicated problems, such as the problem (1.11), the basic algebraic relations can 
be solved in terms of a generalized Fourier integral. This technique in generic, in the sense 
that it can also be used for the solution of simple problems. 

We use such a simple problem, namely the symmetric Dirichlet problem, to illustrate this 
approach: It is shown in Section 4 that the integral defining ^{ak) can be solved for qN{s). 
For A > 0, qN{s) is given by 

Replacing in this equation "^(dk) with the expression obtained by solving equation (1.17) 
for \I/(q;A;), it follows that q'Ar(s) involves a known integral, as well as an integral containing 
the unknown function "^{k). However, using the analyticity properties of the integrant of 
the latter integral, it can be shown that this integral can be computed in terms of residues. 
Furthermore these residues can be explicitly calculated in terms of the known function G{k). 

The situation for more complicated problems is similar: The unknown boundary val- 
ues can be expressed in terms of known integrals, as well as integrals containing the three 
unknown functions {'$.j{k)}l. Exploiting the analyticity properties of the integrants of the 
latter integrals, it can be shown that these integrals can be computed explicitly. 



1.4 Integral Representations for q{x,y) 

When both the Dirichlet and the Neumann boundary values are known, the solution q{x, y) 
can be determined either using the classical integral representation in terms of Green's 
functions [3], or using the novel integral representations constructed in [5] and [8]. For 
completeness, both representations are presented in Section 5. 
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1.5 Organization of the Paper 



In Section 2 we derive the global relation (1.13). In Section 3 we solve the symmetric 
Dirichlet problem (Proposition 3.1), the general Dirichlet problem (Proposition 3.2), the 
general Neumann problem (Proposition 3.3), and we also discuss the oblique Robin problem. 
In Section 4 we discuss the basic algebraic relations associated with the Poincare boundary 
condition (1.11) and derive the relations (1.12b) and (1.12c). In Section 5 we obtain an 
alternative representation for the symmetric Dirichlet problem and then analyze the problem 
defined by equations (1.11) and (1-12). A particular case of this problem, which is solved 
in detail in Proposition 5.1, is a mixed boundary value problem for the modified Helmholtz 
equation. In Section 6 wc discuss the associated integral representations for q[x,y). Further 
discussion of these results is presented in Section 7. 



2 The Global Relation 

Writing the basic elliptic equation (1.7) in the complex variables {z, z) we find 

q,, - Ag = 0. (2.1) 
It is straightforward to verify that this equation can be rewritten in the form [5] 

exp < —ikz — T^l^^) + ^ ( cxp < —ikz — q] — 0, (2.2) 



ik I J - ik \ I ik 

where for the rest of this section k E C — {0}. Suppose that equation (2.1) is valid in a 
simply connected bounded domain Q G C with a piecewise smooth boundary dfl. Then 
equation (2.2) and the complex form of Green's theorem imply 

f exp \ —ikz -^1 ( qzdz rqdz] — 0. (2.3) 

Jan I «^ J V ik J 

In the particular case that Q is the triangular domain D, equation (2.3) becomes 

3 

5^p,(fc) = 0, (2.4) 

where the function Pj{k) is given by the following line integral along the side (j) of the 
equilateral triangle 

Pj{k) = exp !^-ikz - ^z^ (^q^dz - jj^Qdz^ , j = 1, 2, 3. (2.5) 

In what follows we will show that 

pi{k) = pi{k), P2{k) = p2{ak), p-i{k) = p^{ak), (2.6) 
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where the functions Pj{k) are defined in terms of the functions ^j{k) and "ifjik) by the 
equation 



pj{k) = E{-ik) 



, J = 1,2,3. 



(2.7) 



For this purpose we will use the following local parameterizations: 
Side 1: On the side (1) the variable z can be parameterized as 



z{s) 



I 



2^3 



-\-is, s e 



l_ l_ 

2' 2 



(2.8a) 



Then z{—l/2) = Z2 and z{l/2) — Zi. Since the normal and the tangential derivatives are 
parallel to the x and to y axes respectively, it follows that 



^^ = \ i^x - idy) = ^ (9jv - idr) ■ 



(2.86) 



Side 2: If z varies along the side (2) and C varies along the side (1), then z — Cexp {—i^^- 
Thus 



z{s) 



+ is \ e ^ , s e 



2^3 



2' 2 



(2.9a) 



Note again that z{—l/2) — Z3 and z{l/2) — Z2. The equation dz — exp {i^} 5^ imphes that 

(2.96) 



a 



Side 3: In analogy with equations (2.9), if z varies along side (3) we find the equations 



z{s) 



I 



2^3 



and 



zs e 3 , s e 



a 



L L 
'2' 2 



= ^ {dN - idr) ■ 



(2.10a) 



(2.106) 



Finally, z{—l/2) = zi and z{l/2) — Z3. Using equations (2.8)-(2.10) in the expressions (2.5) 
we find equations (2.6), (2.7). 

3 The Analysis of the Global Relation for Simple 
Boundary Value Problems 

3.1 The Symmetric Dirichlet Problem 

We first give the details for the symmetric problem. In this case 



q'^^\s) = qN{s), ^j{k)=F{k), ^j{k) = ^{k), j = 1,2,3, 



(3.1) 
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where the function ^(/c) is defined in terms of the unknown function qN{s) by equation 
(1.14a) (without the superscript (j)), and the function F{k) is defined in terms of the given 
boundary condition /(s) by equation (1.14b), i.e., by the equation 



exp ^ [ k + ^] s 



ds, keC- {0}. 



(3.2) 



Using equations (3.1), the global relation (1.13) and its Schwarz conjugate yield (1-15) and 
(1.16), with 

A{k) = e{ak)F{k) + e{-k)F{ak) + F{ak), 
B{k) = e{-ak)F{k) + F{ak) + e{k)F{ak). 
Hence, since G — A + B, 



G{sn) — 2 cosh 



aSn + — 



A \ / 



aSn I 2 



F{sn) + 2e^""F(as„) + 2F(as„). 



(3.3) 



In summary, we have derived the following result: 



Proposition 3.1 Let the real valued function q{x^ y) satisfy equation (1.7) in the triangular 
domain D, with the Dirichlet conditions (1.8), where 



fj{s)^f{s), 1,2,3, se 



I I 
2' 2 



(3.4) 



and the function f{s) is sufficiently smooth and satisfies the continuity condition f{—l/2) = 
/(//2). Then the Neumann boundary values are the same, q^^is) = qN{s), j = 1,2,3, and 
are given by the Fourier series 



qN[s 



oo 

t X ^ 2in7TS 

r ' 

— oo 



) = 7E 



G(s„) 



sinh 



aSn + 



A \ I 



(3.5) 



where s„ is defined by equation (1.18) and G(sn) is given in terms of f{s) by equations (3.2) 
and (3.3). 



3.2 The General Dirichlet Problem 

The global relation and its Schwarz conjugate yield equations (1-15) and (1.16), where the 
known functions A{k) and B{k) are now given by the equations 



where 



A{k) = e{ak)Fi{k) + e{-k)F2{ak) + F^{ak), 
B{k) = e(-aA;)Fi(A;) + F2{ak) + e{k)F^{6ik), 



(3.6) 



exp { ( ^ + ^ ) 



ds, j = l,2,3, A;eC-{0}. (3.7) 
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Replacing in equations (1-15) and (1-16) k by ak and then eliminating 'i>i{ak) from the 
resulting two equations we find 

e{-ak)'^3{k) + e{k)'^2ic(k) - 2ie{-ak)A{ak) 

= e{ak)'^2{k) + e{-k)'^3{o(k) + 2ie{ak)B{ak). (3.8) 

Taking the Schwarz conjugate of this equation (or cquivalently eliminating ^i(aA;) from the 
equations obtained from equations (1-15) and (1-16) by replacing k with ak) we find 



e{-ak)'^3{k) + e(A;)^2(aA;) + 2ie{-ak)A{ak) 



= e(ak)'^2(k) + e(-k)'^3{ak) - 2ie(ak)B(ak). (3.9) 

Substituting ^2(0;/?) from equation (3.8) and ^3(0;/;;) from equation (3.9) into equation (1.16), 
we find an equation involving \&3(q;/c), ^2(<5A;), and {'^j{k)}i. Eliminating ^3(0;^) from this 
equation and from equation (1.15) we find equation (1.20) with X{k) given by the following 
equation, 



(3.10a) 



X(k) = [e'^(-k)e(ak) + e(-ak)] [Fi(k) + e'^(k)Fs(k)] 
+e\-k) [e^{-k)e{ak)e^{k) + e{-ak)] F^ik) 
+2e'^(k)Fi(ak) + 2F2{ak) + 2e'^(-k)Fs(ak) 

+e^{-k) [2e'^{k)Fi{ak) + {e%k) + l)F2{ak) + 2e^{-k)e\k)F^{ak)] . 
Letting k — kn we find 

X{kn) = [e^{-kn)e{akn) + e{-akn)] x 

[Fiikn) + e\-kn)F2ikn) + e\kn)F^{kn)] 

+2e\kn) [Fi{akn) + e\-kn)F2{akn) + e\kn)F3{akn)] 

+2e{-kn) [Fi{akn) + e^{-kn)F2{akn) + e^{kn)F^{akn)] 
Solving the algebraic equations (1.23) we find the following result: 

Proposition 3.2 Let the real valued function q[x, y) satisfy equation (1.7) in the triangular 
domain D, with the boundary conditions (1.8), where the given functions /j(s) have sufficient 
smoothness and are continuous at the vertices. Then the Neumann data q^^{s)., j = 1,2,3 
can be expressed in terms of the given Dirichlet data by the Fourier series 

CXD 

<l^\s) = ^ E [^(^3n) + C?)e^M(^3n-l) + ^'^e^^ M(^3„-2)] e"^, (3.11) 



(3.106) 
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where km is defined by equation (1.21), 



.(1) _ Ji) 



2 — 

M{km) = 



(2) _ ^(3) 



--2 — '-^1 "^1 
2iX{km) 



.(3) _ (2) 



a. 



a'^e(akm) - e(-akm) 
and X{km) is defined in terms of fj{s) by equations (3.7) and (3.10). 

3.3 The Oblique Robin Problem 

Suppose that q{x,y) satisfies the Poincare boundary condition (1.11), i.e. 



(3.12) 
(3.13) 



(3.14) 



sm Pj \ " ds 

Substituting this expression in the definition of Pj{k), i.e. in the equation (2.7), we obtain 



Pj{k) = E{-ik) I ^ cxp <! l^/c + ^ I s 



1 d 



A 



ds. 



Integrating by parts we find the following expression for Pj{k): 

p,{k) = iE{-ik) [H,{k)Y,{k) + F,{k) + C,{k)] , j = l, 2, 3, 
where the function Hj{k) is defined by 

A 



Hj{k) = ke'^^ + 



(3.15) 



(3.16) 



the function Fj{k) is defined in terms of the given boundary conditions fj{s) by the equation 



2 sin 13 j 



exp <j [k + ^]s\ fj{s)ds, 



(3.17) 



the function Cj{k) involves the values of q{x,y) at the vertices. 



e{—k)q 



e{k)q^^^ ( ^ 



2 sin Pj 

and the function Yj{k) involves the unknown Dirichlet boundary values. 



(3.18) 



Y^k) 



1 



2 sin 

In equations (3.15)-(3.19), k is complex and k 0. 



cxp ^ [ k + j^] s\ q^^\s)ds. 



(3.19) 
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In the particular case of the obhque Robin problem (1.9), (5j — (5 and 7^ =7, j — 1, 2, 3, 
thus Hj{k) = H{k), where H{k) is defined by equation (3.16) without the subscript (j). 
Substituting the expression for Pj{k) (with Hj — H) in the global relation (2.4) we find 

E{-ik) [H{k)Yi{k) + Fi{k) + Ci(A;)] + E(-iak) [H(ak)Y2(ak) + F2(ak) + C2(ak)] 

+E{-iak) [H{ak)Y^{ak) + F^{ak) + C^{ak)] = 0. (3.20) 

The contribution from the corner terms Cj cancels. Indeed, this contribution is proportional 
to the following expression 



E{-ik) 



e{-k)q 



(1) 



e(^),(^) ( L 



+E{-iak) 



e(—ak)q 



(2) 



e{ak)q^'^ ( ^ 



+E{-iak) 



-ak)q 



(3) 



e{ak)q^'^ ( ^ 



But, the assumption of continuity at the vertices implies 



,(1) 



r(2) 



.(2) 



,(3) 



,(3) 



,(1) 



Hence the terms q'^^\—l/2) and q^'^\l/2) in the expression (3.21) cancel iff 

E{-ik)e{-k) = E{-iak)e{ak). 
This equation is indeed valid, and it is the consequence of the identity 



(3.21) 



(3.22) 



(3.23) 



1 

7! 



-ik + 



-ik 



+ -k + 



1 

7! 



-iak + 



X 



-iak 



+ { ak + 



ak 



Using the fact that the corners term cancel, the global relation and its Schwarz conjugate 
yield (compare with equations (1.15) and (1.16)) the following equations: 

e(ak)H(k)Yi(k) + e{-k)H{ak)Y2{ak) + H{ak)Y^{ak) = -A{k), 

e{-ak)H{k)Yi{k) + H{ak)Y2{ak) + e{k)H{ak)Ys{ak) = -B{k), (3.24) 

where A{k) and B{k) are defined by equations (3.6) in terms of Fj. 
Let 

P{k) = £ffi. (3.25) 
H{k) 

Following precisely the same steps used for the general Dirichlet problem we find the following 
expression for Y2{ak) in terms of {Yj{k)}l: 



,.P\ak) 3,,,P{ak) 



P^{ak) 



P{ak) 



H{ak) 



Y2{ak) = T{k) 



+ 



■ , 2, .^p{ak)p{k) , 

e{-ak) - e\-k)^^^l-^e{ak) 
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Pjak) 
P{ak) 



e{—ak) — e^{k) 



P^{ak) P{ak)P{k) 



P^{ak) P\ak) 
where T{k) is given in terms of Fj{k) by the following equation 

P{ak) 



e{ak) 



Y2{k), 



(3.26) 



+ 



H{k)T{k) = e{-k) 



E\-iak)^^- E\iak) ^ 



E\-iak) - E\iak) 



P^{ak) 



P{ak) 



P{ak) 



F2{k) + e{k) 



E\-iak)^^- E\iak) ^ 



P{ak) 



+e\k)^—^-—Fi{ak)+ ^ ^ 



P{ak) 



P{ak) 



P^{ak) 



P{ak) 
Fsiak) 



+e{-k) 



^ ' Fi{ak)+[e^{k) ^ ' 



P{ak) ' V P{ak) 
li k — kjn, where km is defined by 



e (-^)l^7rVv F2{ak)+e{k) Fs{ak). 



P^{ak) 



P{ak) 



(3.27) 



' P^iak) ~ 



then 



2i7rm . . , , 

= 63, e {-k) 



2, j^P{ak) 2inm 



— e 3 



P{ak) ' ' 'P{ak) 

Thus evaluating equation (3.26) &X, k — km, we find the following expression: 



where 



G{k, 



2T(fc,„)sin/? 



G{km), meZ (3.28) 



(3.29) 



Even if one is able to invert equation (3.28) for the bracket appearing in the integrant 
of the Ihs of equation (3.28), one will find an expression which will involve an infinite series 
over the transcendental values of km- Thus instead of analyzing the relevant inversion we 
will analyze the general Robin problem using the generalized Fourier transform approach 
(see Section 5). 

In the case of the Neumann problem, equations (3.28) and (3.29) simphfy and yield the 
following result. 



Proposition 3.3 Let the real valued function q{x, y) satisfy (1.7) in the triangular domain 
D, with the Neumann boundary conditions 



/ / 

2' 2 



J = l,2,3 



(3.30) 
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where the functions fj{s) have sufficient smoothness and are continuous at the vertices of 
the triangle. Then the Dirichlet data q^^\s), j — 1,2,3 can be expressed in terms of the 
given Neumann data by the Fourier series 



oo 

3i i: 



^•^ ' £3 'XI 



e ' 



(3.31) 



where c^p are given by (3.12) and 



a'^e{akjn) — e{—akm) 
The known function Tr^{km) is defined by the equation 
A 



(3.32) 



ik + — ) TN{k) = e{-k) [E^{-iak) + E^{iak)\ Fi{k) + [E^-iak) + E^{iak)] F2{k) 

+e(k) [E^(-iak) + E^(iak)] F^ik) + 2e^{k)F^{ak) + 2F2{ak) + 2e'^{-k)F^{ak) 

+2e(-k)Fi(ak) + {e^k) + e^(-k)) F2{ak) + 2e{k)F^{ak), (3.33) 
where Fj{k) is given by (3.17) with sin/?^ = 1. 



4 Poincare Type Boundary Value Problems 

Suppose that q{x, y) satisfies the Poincare type boundary condition (1.11). Then substituting 
the expression Pj{k) from equation (3.15) into the global relation (2.4), we find an equation 
similar with equation (3.20), where H{k), H{ak) and H{ak) are replaced by Hi{k), H2{ak) 
and Hslak), and Fj, Cj, Yj are defined by equations (3.17)-(3.19). Proceeding as in Section 
3.3, in analogy with equation (3.26), we now find 



D{k)H2{ak)Y2{ak) = ^Tj{k)Hj{k)Yj{k)+T{k)^C{k), 



(4.1) 



where T{k) is defined in terms of the known functions /j(s), C{k) involves the values of q at 
the corners, and D{k), {Tj{k)}\ are defined by the following equations: 



Pi{ak) 



V,{k) 



P2{ak)P-i{ak) 
1 



3^ 3n.^Pi{cik)P2{ak)P3{ak) 



e\-k) - e\k) 



V2{k) = e\-k) 

r,{k) = e\k) 



Pi{ak) 



e{—ak) — e {—k)e{ak) 



Pi{ak)P2{ak)P3{ak) 
Pi{k)Pi{ak) 



P2{k)P2{ak) 
Pi{ak) 



P3{k)Ps{ak) 



e{—ak) — e^{k)e{ak) 
e{—ak) — e'^{—k)e{ak) 



P2{ak)P3{ak)_ 

P2{k)P2{ak) 



Pi{ak)P3{ak) 

P3{k)P3{ak) 
Pi{ak)P2{ak) 



(4.2) 

(4.3a) 
(4.36) 
(4.3c) 
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with 

Pj{k) = SM. (4.4) 

In order to be able to solve this problem using a generalized Fourier integral we require 
that when D{k) vanishes, then T2{k) and r3(A;) are proportional to ri(A;). Actually, T^i^k) 
is proportional to ri(A;) for all complex k provided that 

Pi{k)Pi{ak)Pi{ak) = P^{k)Pi{ak)P^{ak). (4.5a) 

Equating the brackets appearing in the definitions of ^i{k) and T2{k), and replacing in the 
resulting expression e^{k) by 

Pi{ak)P2{ak)P:i{ak) 
Pi{ak)P2{ak)P3{ak)' 
it follows that r2(/c) is proportional to ri(fc) provided that 

Pi{k)Pi{ak)Pi{ak) - P2{k)P2{ak)P2{ak). (4.56) 

Equation (4.5b) is valid if the following two equations are valid: 

sin3(/5i-/32) = 0, 72(3A-7|)sin3/?i-7i(3A-72) sin 3/32 = 0. (4.6) 

Indeed, in order to simplify equation (4.5b) we first compute the product 

Hi{k)Hi{ak)Hi{ak), 

H,{k)H,{ak)H,{ak) = k^e''^' + + 3A71 - 7?. (4.7) 



The function Hi(k) can be obtained from Hi{k) by replacing f3i with — /3i, thus Hi(k) Hi(ak) 
Hi(ak) is given by an expression similar to (4.7) with (3i replacing by Hence equation 
(4.5b) yields 

f + kJ^TiTlT + 'J^ ;1 ^ 7i A f - + + 'J^ ;2 - I2 



k'e-^'^^ + + 3A71 - 7? k^e-^'^^ + + 3A72 - 7f ' 

This equation simplifies to the equation 

(^k^ - sin 3(/3i -P2)+(^k'- [(3A72 - 72') sin 3/3i - (3A71 - 7^^) sin 3/^2] = 0, 

which is valid for all k iff equations (4.6) are valid. 

Equation (4.6a) implies P2 = Pi + nn/3, then sin 3/^2 = sin 3/9i exp{m7r} and we find 
equation (1.12b). Similarly equation (4.5a) yields equation (1.12c). 

The Case that the Corner Terms Cancel 

The definition of the corner terms Cj{k), i.e. equation (3.18), shows that Cj{k) involves 
exp[iPj]/ sinPj. Thus the contribution of the corner terms in the global relation (3.20) 
vanishes iff 
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Example 1. 



In this case 



A = /32 = /33 = arbitrary. 



ka 

Example 2. 

A = /92 = /93 = A /3 arbitrary, 72 = 73 = 0, 71 = (3A)^ A > 0. 
In this case 



In particular if /3 = 7r/2, then 



A 



ke'l^' 



H,{k) = i ( A; - - ) - 71, H2{k) = Hs{k) ^i{k- 



Thus 



Hence 



~^ \.'^~ k) ~ 



1 



r2(A;) = e^{-k)Pi{ak) 



Pi{ak)_ 

[e{-ak) - e^{-k)e{ak)Pi{k)Pi{ak)] , 

e'^{k)e(dk) 



r^ik) = e^{k)Pi{ak) 



e{—ak) + 
e(— a/c) + 



e2(-A';)e(«A;) 



Pi{ak) 

If -Pi(A;) is defined by (4.14) with 7^ = 3A, it can be verified that 

1 



Pi{k)Pi{ak) 



' Pi{ak)' 



Hence 



e\k) e\-k)T,{k) 
-Vi[k), T2[k)\D{k)=Q = - 



Pi{aky 



Pi{ak) 



The Laplace Equation 

In the particular case of the Laplace equation with 7^ = 0, it follows that Pj 
Pj is independent of k. 
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5 The Analysis of the Global Relation Via Fourier In- 
tegrals 

In this section we restrict A to be non-negative. Slightly more complicated formulae can be 
derived for A < 0. 

We first derive equation (1.25). Letting k — \k\e^ , the definition of ^{ak) yields 

*(-i|/e|) = f\xp\(-i\k\ + ^--)s\qN{s)ds. (5.1) 



2 



-i\k\ 



Suppose that A > 0. Letting t{\k\) — \k\ — X/\k\, it follows that if \k\ e (0, +00), then 
t e (—00,00). Thus inverting equation (5.1) we find 



1 



qN{s) = — / e'''-^{-i\k\)dt, s e 



L L 
'2' 2 



, ?jv(s) — elsewhere, 
where (in the argument of ^) is a function of t. Rewriting t in terms of \k\ we find 



Letting \k\ — ke ''^/^ we obtain 



Qn{s) — — / exp |— ( ak + — ] si ^(a/c) ( ke' 
2n Jq I \ akj } \ 



X \ dk . ^, 



The rhs of this equation equals 



27r J_(3oe^ 



exp ^ - ( + ^) ^(aO (^e"^ + ^) |. (5.3) 



Indeed, for the derivation of (5.3) we first observe that the function '^{ak) remains invariant 
under the transformation k — > aX/k. Thus making the change of variables k — aX/ ^ in the 
rhs of equation (5.2), and using 

, _HL A /, _i7L X \ dk d$ 

« + ^ ^ - Ue ^ + 



we find the expression (5.3). Combining (5.2) and (5.3) we obtain 



1 ( ( , A \ 1 _ , . /. X \dk 



Using — ia, this equations becomes equation (1.25). 
If A = 0, we set k — texp{i^}, t e R, and rewrite (1.14a) as 

^{it) = / e'^'qN{s)ds (5.5) 

J-l/2 
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which is inverted to 

7 r+oo 

QNis) = — / e-'''^iit)dt. (5.6) 
Replacing in (5.6) {it} with {ak}, we arrive at 

= ^ / e-'^'=^*(aA;)ciA;. 

This equation, in comparison to (1.25) misses a factor of 1/2; this is due to the hnearity of 
the relevant transformation in this case. 

5.1 The Symmetric Dirichlet Problem 

Solving equation (1.17) for "ifi^ak) and substituting the resulting expression in equation (1.25) 
we find 

(5.7) 

where A(A;) denotes the coefficient of ^{ak) in equation (1.17), i.e. 

A{k) = e{k)-e{-k). 

The line ^ooe"^, ooe^^ splits the complex /c-plane into the two half planes, 

V+^heC, ^<argA;<^|, 



^ I , ^ 77r , IStt 

V-^{keC, — <arg/c< — 



We observe that 



exp |— ^ak + — ^ s I e(aA;) is bounded for k eV , 

exp I— ^a/c + — ^ e(— a/c) is bounded for k&V~^. (5.8) 

Indeed the exponential of (5.8a) involves ak (| — -s), and since 1/2 — s >Q, the exponential 
in (5.8a) is bounded for Re{ak) < 0, i.e. in V". Similarly, the exponential of (5.8b) involves 
—ak (I + s), which since //2 + s > 0, is bounded for lie(ak) > 0, i.e. in T)^. 

We also note that '${k)/A{k) is bounded for all k ^ C, k ^ s„. Indeed, for Re /c > 0, 
A(A;) is dominated by e{k), while for Re /c < 0, A{k) is dominated by e{—k), hence 



^{k)e{-k), Rek > 

(5.9) 

-*(A;)e(A;), Rek < 0. 
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Furthermore ^(/c)e(— A;) involves k{s — 1/2) which is bounded for Re A; > 0, while ^{k)e{k) 
involves k{s + 1/2) which is bounded for Re A; < (recall that — | < s < |). 

The above considerations imply that the parts of the integral (5.7) containing e{ak)'^{k) 
and e(— q:A;)^(A;) can be computed by using Cauchy's theorem in V~ and respectively. 
The associated residues can be computed as follows: Let s+ and s~ denote the subsets of s„ 
in D"*" and T>~, respectively. Evaluating equation (1.17) at A; = we find 



2iG{s; 



-e^{-as-) + 1' 



-e(-a<)*(4) 



2iG{s; 



1 — e^^Q;s^ 



Thus 



—ia exp < — 



- { ak + — 
ak 



G{k) 

m 



A 



A'{s+) [l-e2(asH 



(aA;)2 
1 - ^ 



dk 



A 



[as2 



+ia exp < — ( aSj^ + 



as: 



Gis-) 



1 - 



A 



(5.10) 



5.2 The Poincare Problem 

Evaluating equation (4.1) at A; = km, where km is a zero of D{k), it follows that the unknown 
terms Yj{k) appear in the form 

2/> ^ Pi{akm)Pi{km) H^{km) ^ f,^ 

Pi{km)Pz{0lkm) Hi{km) 

The crucial difference of this general compared with the oblique Robin case (1.9), 

is the following: Using the definition of Yj{km) we find that the coefficients of q^'^\s) and of 
q^^\s) involve in general Aim-dependent expressions, thus it is not clear how the associated 
integral can be inverted. In contrast, equation (4.1) can be solved using the approach of 
section (5.1). The definition of y2(Q;A;), i.e. equation (3.19), and equation(1.25), imply 

w = ™ £j {- - s) 4 - wp) 

Solving equation (4.1) for Y2{ak) and substituting the resulting expression in equation (5.11) 
we find an integral involving the three unknown functions {Yj{k)}l. The unknown part of 
this integral involves the factors (5.8) analyzed already, as well as factors of the type 

^ ^D{ky D{k) ■ 
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These terms are bounded for all k ^ km- Indeed, ignoring the terms involving Pj{k) we find 

Yj{k)e{-k) -e^-k), Rek > 

Yj{k)e{k), Rek < 0, 

Yj{k)e{-k), Rek > 



e\k) 



D{k) 



Yj{k)e{k)-e^{k), ReA; < 0, 

which are identical with the expressions (5.9) except for the occurrence of the factors e^{—k) 
and e*^(A;) for Re /c > and Re /c < 0, which are bounded. 
The above discussion implies that the integral involving 



-ak) 



Hijk) 



Yi{k) 



-k) 



P2{k)P2{ak) 



Ps{k)Ps{ak) 



Ys{k) 



(5.12) 



can be computed by using Cauchy's theorem in 1?+. Evaluating equation (4.1) at /c+ it 
follows that the associated residue equals 



-[T(A;+)+C(A;+)] 



(5.13) 



P2{ak+)P3(,ak+) 



Similarly, the integral involving 

Pi{k)Pi{ak) e{ak) 



P2{ak)P^{ak) D{k) 



-k) 



H,{k)Y,{k) + 



e^{k)P2{ak)P^{ak) 
Pi(k)Pi(ak)P3(ak) 



H2{k)Y2{k) 



Pz{ak) 



:Hs{k)Ys{k) 



Pi{ak)Pi{ak)"^^'"'^''^'"'\ ^^'^^^ 

can be computed by using Cauchy's theorem in . Evaluating equation (4.1) at /c^, it 
follows that the associated residue equals 

\nk-j + ciK 



In what follows we give the details for a mixed Neumann-Robin problem. 



(5.15) 



We will consider Example 2, as it is described by (4.11) with /3 = |. On side (1) we 
assume the Robin condition 

(5.16) 



q^i,\s) + V^q^'\s) = h{s), 
and on sides (2) and (3) we assume the Neumann conditions 

<lf{s)=f2{s) 

and 

<lf{s) = h{s). 



(5.17) 
(5.18) 
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Then Hj{k), Pj{k), D{k) and Vj{k), are given by (4.13), (4.14), (4.15a) and (4.15b,c,d) 
respectively. Furthermore, T^{k) is proportional to Ti{k) (see equation (4.16a)), while r2(A;) 
becomes proportional to ri(A;) only on those /c^'s for which D{k) vanishes. These are roots 
of the transcendental equation 



exp 



(fc+|-%A)(A: + A-2VA) 



V\]{k + ^ + 2V\ 



(5.19) 



For k — km, equation (4.1), in view of (4.16), implies 



Tl{km) 



Hl{km)Y,{km) 



-H2{km)Y2{km) 



Pi{akm) 



Pi{akm) 

By virtue of (4.15b) and the identity 

Pi{km)Pi{akm)Pi{akm) = -1 

equation (5.20) is written as 



Hs{km)Ys{km) 



-T{km). 

(5.20) 
(5.21) 



ei—akjn) + eiakn 



C^i kr, 



Pl{km) 



Pi{akm)_ 

+e^{-k„,)Pi{ak^)H2{kjn)Y2{km) + e^{k^)Pi{akjn)H^{km)Y'i{k„^)] = -T{k„,). (5.22) 
Since the corner terms C{k) vanish, the representation (5.11) and equation (4.1) yield 



2^ iooe*^ ^""^ 



a/c + ^ 1 s H 1 



A 



1 



J2^j{k)Hj{k)Y,{k)+T{k) 



H2{ak)D{k) 

Utihzing the expression (5.22) we arrive at the following result. 



{aky 



dk. 



(5.23) 



Proposition 5.1. Let the real valued function q{x,y) satisfy equation (1.7) with A > 
in the triangular domain D, with the Robin boundary condition (5.16) on side (1) and 
the Neumann boundary conditions (5.17) and (5.18) on sides (2) and (3), where the given 
functions fj{s) have sufficient smoothness and are continuous at the vertices. Then the 
Dirichlet value on side (2) is given by 



ooe D 



exp 



ak + — I s 

ak 



A 



k+ 



exp 



H2{akt,)D'{kt,) 



A 



aky J H2{ak)D{k) 



dk 



Pi(afc^) 
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exp{-(av;.H-^),} . ^ 

^ H2(a*^-)£>'(A-) V (ak-)^ J 1 + Piiak-JE'i-mk-y ^' ^ 



^77 

where D'{k^) denotes the derivative of D{k) evaluated ai k = k^. The summations are 
taken over all k^ G V'^ and A;~ G 7?" respectively, and T, H2, Pi are defined by equations 
(3.27), (3.16), (4.4) respectively. 

There exist similar formulas for the Dirichlet values on sides (1) and (3). 

6 The Integral Representations 

If A > the classical Green's representation is given by [3] 

g(r) = ^ / \K{2VX\r - r'|)a„,g(r') - q{r')dn>K{2VX\r - r'|)l dZ(r') (6.1) 

^TT Jqd L J 

where the integration is over the boundary dD of the triangle in the positive direction, dn' 
denotes the outward normal derivative on dD, dl{r') is the line element along dD, and K{x) 
is the modified Bessel function of the zeroth order and of the second kind for the modified 
Helmholtz equation. For the case of the Helmholtz equation K{x) is proportional to the 
Hankel function of the zeroth order and of the first kind, while for the Laplace's equation 
K{x) is proportional to the logarithm of x. 

For the Laplace equation, the integral representation constructed in [5] is defined as 
follows: ^ 

|! = ^E / ^'''Pjik)dk, zeD, (6.2) 
j=i -^ij 

where the contours Ij are the rays from to 00 specified by the arguments — 7r/2, 7r/6, 57r/6 
respectively, and the functions Pj{k) are defined by equations (2.6) in terms of Pj{k), where 
the latter functions are defined by equations (2.7), (1.14) with A = 0. 

ki 




Figure 6.1: The rays Ij in the complex /c-plane 
For the modified Helmholtz equation, the analogue of equation (6.2) is [5] 



q[z 



e'''+r,--pj{k)'^^ 



k ' 



zeD, 



(6.3) 
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where the rays Ij are the same as in (6.2) and Pj{k) are defined by equations (2.6) and (2.7). 

There exists a similar representation for the Helmholtz equation, which however, in ad- 
dition to rays, it also involves circular arcs [5]. 

6.1 The Symmetric Dirichlet Problem 

Using the integral representation (6.3) it is possible to compute directly q{z,z), bypassing 
the computation of the unknown boundary values. For brevity of presentation we will only 
give details for the symmetric Dirichlet problem. The analysis of the more general boundary 
value problems (1.8)-(1.11) is similar. 

Recalhng the definitions of pj, i.e. equations (2.6) and (2.7), it follows that the represen- 
tation q{z, z) given by equation (6.3) involves the known function F{k) defined by equation 
(3.2), as well as the unknown function ^{k) which on the rays Ij appears as: 

h ■■ E{-ik)'^{k), h : E(-iak)^(ak), h : E{-iak)'^{ak). (6.4) 

Solving equation (1.17) for '^{ak) in terms of ^(A;), and then using the Schwarz conjugation 
of the resulting equations in order to express '^{ak) in terms of ^{k), it follows that the 
expressions in (6.4) involve the unknown function ^(A;)/A(A;), A(A;) = e{k) — e{—k), times 
the following expressions: 

li : E{—ik)[e{k) — e{—k)], I2 ■ E{—iak)[e{ak) — e{—ak)], I3 : E{—iak)[e{ak) — e{—ak)]. 

(6.5) 

The third of the relations in (1.5) implies e{k) — E{iak)E{—iak), thus 
E{-ik)e{k) = E\iak), E{-ik)e{-k) = E^{iak). 
Replacing k by ak and by ak in these identities, we find 

E{-iak)e{ak) = E^{iak), E{-iak)E{-ak) = E^{ik), 

E{-iak)e{ak) = E'^{ik), E{-iak)e{-ak) = E^{iak). 
Thus the expressions in (6.5) involve 

h : E\iak) - E^{iak), h : E'^iiak) - E\ik), h ■ E\ik) - E^{iak). 

Hence, the unknown part of q{z, z) involves the following integral 

3 

J(^^^) ^^Jj(^^^)^ (6-6) 

M^^^)^^f exp likz + ^z\ E^{iak) , (6.7a) 

J,{z,z)^^ [ expUkz + ^z\E\ik)^^^^, (6.76) 
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{-«3}U{/l} 



exp < ikz H — -z > E iiak) , . . . 
1^ z/c J kA{k) 



(6.7c) 



Each of the above integrals can be computed in terms of residues. Indeed, it was shown in 
Section 5 that '^{k)/kA{k) is bounded as k ^ and as /c — > oo. Furthermore, it will be 
verified below that the exponentials. 



A 



ik 



A 



ik 



exp < ikz + —z > E (iak), exp < ikz + —z > E (ik), exp < ikz + —z > E (iak), (6.8) 



A 



ik 



are bounded as /c — > and /c — > oo, for arg k in 



" TT TT" 




TT 57r 




57r 


37r" 


.~2' 6. 




6' Y 






Y 



respectively, provided that {z,z) G D. We first consider the first exponential in (6.8); since 
Z2 — this exponential can be written as exTp{ik{z — Z2) + X{z — Z2)/ik}. If z is in the 

triangular domain then 

TT StT 

- < arg(^ - Z2) < — . 
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Thus if — I < arg A; < |, we find 



< arg[A;(2; — Z2)] < tt. 



Hence exp{ik{z — Z2)} is bounded as \k\ — > 00 and exp{A(^ — Z2)k/i\k\'^} is bounded as 
\k\ — > 0. Hence there is no contribution from zero and from infinity. 

The results for the second and for the third integrals in (6.8) follows from the above result 
by using appropriate rotations. The roots of A (A;) = lie on the imaginary axis. Denote 
by those with positive imaginary part and by those with negative imaginary part. 
Obviously, the residue from each s+ has a full contribution to Ji{z,z), while the residue 
contribution from each s~ is split in two halfs, one half is contributed to J2{z,z) and one 
half to J3{z, z). 

Tedious but straightforward calculations lead to the expression 



E\ist)G{sl) 



s+A'(s+)A(a4) 

E\ias-) + E\ias-))G{s-) 
«nA'(s-)A(as-) 



(6.9) 



In the above calculations, the value of ^i(s^) is obtained from (1.18). 
We observe that for each n e Z, 



A 

exp \ ISnZ + - — z 

%Srt. 



exp |±2y(^)2 + Aa;| • exp {-2^^y} (6.10) 



thus this expression shows that the equilateral triangle admits separable solutions. It is clear 
that each eigensolution in (6.10), solves equation (1.7). 
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7 Conclusion 



Eigenvalues and eigenfunctions for equation (1.7) with homogeneous Dirichlet, Neumann, 
and Robin boundary conditions were constructed in the classical works of Lame [10]- [12]. 
Some of these results have been rederived by several authors, in particular the Dirichlet 
problem is discussed in the recent review [13]. The Robin problem is analysed in [14]. It is 
remarkable that Lame argued, using physical considerations, that it is impossible to solve 
certain problems using infinite series as opposed to integrals. Indeed Lame writes [12 p. 191]: 
"The series should therefore express the fact that the temperature remains zero on strips 
of constant width separated by other strips of double width, in which the temperature may 
vary. The analytic interpretation of this sort of discontinuity demands the introduction of 
terms where the variables appear inside integrals^. These terms, of a nature that we will not 
consider here, cannot disappear from the total series unless the discontinuity disappears" . 

In this paper we have solved several boundary value problems by introducing a novel 
analysis of the global relation, i.e. of equation (1.13). Although this equation was first 
derived in the important work [14], where it was also used to solve the Robin problem, our 
treatment of equation (1.13) is different than that of [14]. As a consequence of our novel 
analysis of equation (1.13) we are able to first present a straightforward treatment of simple 
boundary value problems. This treatment, which is based on the evaluation of the basic 
algebraic relations (see the Introduction) at particular values of k, expresses the unknown 
boundary values in terms of infinite series. The Dirichlet, Neumann and Robin problems can 
be solved using this approach. We then show that, in agreement with the above remarks of 
Lame, more complicated boundary value problems apparently require the use of generalized 
Fourier integrals as opposed to infinite series. Proposition 5.1 presents the solution of such 
a problem. 

In this paper, as opposed to the works of [1], [2], [5]- [7], we have introduced a method 
for determining the generalised Dirichlet to Neumann map, i.e. determining the unknown 
boundary values as opposed to determing q{x,y) itself. In this respect we note that: (a) In 
some applications one requires precisely these unknown boundary values, (b) When both 
the Dirichlet and the Neumann boundary values are known, it is straightforward to compute 
(l{x,y). 

We emphasize however, that the approach of Section 5 can be used to construct directly 
q{x,y). Indeed, if one uses the novel integral representations for q{x,y) obtained in [5], 
instead for the representation (1.25) for g^r, and if one follows the approach of Section 5, one 
can again compute explicitly the contribution of the unknown functions Yj{k). This latter 
approach is illustrated in Section 6.1 for the symmetric Dirichlet problem. More complicated 
problems using this approach are solved in [2], [6], [7]. 

In order to compute q{x,y) from the knowledge of both the Dirichlet and the Neumann 
boundary values one can use either the classical Green's formulae or the representations of 
[5] . Regarding the latter representations we note that they provide a tailor-made transform 
for the particular problem at hand. In fact the exponential exp{ikz + {X/ik)z} reflects the 
structure of the PDE, the contours Ij in the complex fc-plane reflect the geometry of the 
domain, and the functions Pj{k) describe the boundary conditions. 

^ These words were not italic in the original text 
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Both the Dirichlet and the Neumann problems involve elementary trigonometric func- 
tions. It is interesting that the analysis of the global condition yields these separable solutions 
without the direct use of separation of variables. 

For arbitrary values of the constants j3j and 7^, the Poincare problem (1.11) gives rise to 
a matrix Riemann-Hilbert problem. For the particular case that equations (1.12) are valid, 
it is possible to avoid this Riemann-Hilbert problem and to solve the problem in closed form. 
Although equations (1.12) impose severe restrictions on (5j and 7^, some of the resulting 
cases appear interesting. These cases include the following: 

(1) /?! = /?2 = /?3 = 1^, 7j arbitrary 

In this case the angles are specified, but 7^ are arbitrary. 

(2) The mixed Neumann- Robin problem analyzed in Section 5. 

(3) /32=/9 + f, /53 = /5 + f, 71 = 72 = 73 = 7. 

In this case all derivatives are computed along a direction making an angle (3 with the 
positive vertical axis. 

The results presented here can be made rigorous, following a formahsm similar to the 
one used in [9] . 

Several problems remain open which include the following: 

1. The investigation of singularities associated with discontinuous boundary conditions. 

2. If the /5j's differ, then the global relation (3.20) contains a contribution from q at the 
three corners. Several approaches for determining these terms are presented in [1] and [6], 
however, the optimal treatment of these terms remains open. 

The approach introduced in [4] and [5] constructs the solutions of a given boundary 
value problem without using eigenfunction expansions. Similar considerations apply to the 
approach introduced here for constructing the generalised Dirichlet to Neumann map. How- 
ever, it turns out that the above approaches can also be used to investigate the existence of 
eigenfunction expansions and to construct these expansions when they exist. This will be 
presented elsewhere. 
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